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ABSTRACT 

We report on a calculation of large scale anisotropy in the cosmic microwave 
background radiation in the global monopole and texture models for cosmic struc- 
ture formation. We have evolved the six component linear gravitational field along 
with the monopole or texture scalar fields numerically in an expanding universe 
and performed the Sachs- Wolfe integrals directly on the calculated gravitational 
fields. On scales > 7°, we find a Gaussian distribution with an approximately 
scale invariant fluctuation spectrum. The AT/T amplitude is a factor of 4-5 larger 
than the prediction of the standard CDM model with the same Hubble constant 
and density fluctuation normalization. The recently reported COBE-DMR results 
imply that global monopole and texture models require high bias factors or a large 
Hubble constant in contrast to standard CDM which requires very low Ho and 
bias values. For Hq = TO^Mpc" 1 , we find that normalizing to the COBE results 
implies bg ~ 3.2 ± 1.4 (95% c.l.). If we restrict ourselves to the range of bias 
factors thought to be reasonable before the announcement of the COBE results, 
1.5 ^ &8 ~ 2.5, then it is fair to conclude that global monopoles and textures are 
consistent with the COBE results and are a better fit than Standard CDM. 
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Topological defects which formed in a phase transition in the early universe 
provide an attractive mechanism for the generation of density perturbations which 
can grow to form galaxies and large scale structure. Cosmic strings (Vilenkin, 
1980, Zel'dovich, 1980), global monopoles (Bariola and Vilenkin, 1989, Bennett 
and Rhie, 1990), and global textures (Turok, 1989) have all been proposed as 
possible seeds for cosmic structure formation. These theories are characterized 
by a single adjustable parameter, the Grand Unified theory symmetry breaking 
scale, v, and the value (v ~ 10 16 GeV, Gv 2 ~ 10~ 6 ) predicted by particle physics 
(Amaldi et al, 1991) also gives the correct amplitude to generate galaxies and 
large scale structure. In contrast, inflationary models generally cannot produce 
a reasonable perturbation amplitude without a rather extreme fine tuning of the 
coupling constant (A ~ 10~ 12 ). 

Recent studies (Cen et al, 1991, Park et ai, 1991, Gooding et al, 1991 and 
1992, and Turok and Spergel, 1990) of texture seeded density perturbations in a 
universe dominated by cold dark matter (CDM) have indicated that this theory 
may solve two of the perceived problems with the standard CDM model: the 
lack of sufficient large scale structure and quasars at high redshift. Their results 
generally agree with our calculations of global monopole (GM) and texture (T) 
seeded structure formation (Bennett, Rhie, and Weinberg, 1992), but we also find 
that these nongaussian seeds tend to generate large galaxy pair velocities and 
cluster velocity dispersions (see also Bartlett, Gooding and Spergel, 1992). This 
can be alleviated by selecting a larger "bias" factor (i.e., a lower normalization of 
the density field). 

Another major difference between the nongaussian GM & T models and the 
gaussian inflationary models can be seen in the cosmic microwave background 
radiation (CMB) anisotropies. In inflation scenarios, AT/T is due to remnant 
quantum fluctuations crossing the horizon at last scattering, and analytic results 
for AT/T have been obtained for both the scalar mode responsible for the growth of 
cosmic structure (Bond and Efstathiou, 1987) and the tensor modes (gravity waves) 
(Abbott and Wise, 1984). For the standard exponential inflation models, only the 
scalar growing mode is important. With topological defects, metric fluctuations are 
generated by the relativistic dynamics of the defects inside the horizon, and AT/T 
is affected by the fluctuations in all six independent components of the gravitational 
field. Thus, we expect the topological seed models to predict larger AT/T for a 
fixed amplitude of density fluctuations. The recent detection of CMB anisotropy by 
the COBE DMR experiment (Smoot, et ai, 1992) at a level somewhat higher than 
the prediction of the Standard CDM model should be regarded as encouraging for 
topological seed models. (Some non-standard inflationary models may also have 
significant tensor mode perturbations (Davis, et al, 1992), but they are certainly 
not the only theory with tensor modes.) 
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A previous estimation of AT/T for the texture model (Turok and Spergel, 1990) 
was based on a simple analytic model of a single texture evolution. They found a 
non-Gaussian distribution of hot and cold spots at a level that seems to conflict 
with the COBE data. In this paper, we present realistic numerical calculations of 
AT/T on COBE scales in the GM&T models with no simplifying assumptions. 

AT/T on COBE scales reflects the variations in time delay (frequency shift) 
along the photon paths from last scattering until the present. This is the general- 
ized Sachs- Wolfe effect (Sachs and Wolfe, 1967) where not only scalar growing mode 
but all gravitational field components contribute to the temperature fluctuations. If 
we choose a coordinate system such that the metric is g^vix) = a(r]) 2 [r]^ + h^ u (x)], 
where 7]^ v = diag (—,+,+, +) and is the metric perturbation, and choose a 
gauge Hqv = (Veeraraghavan and Stebbins, 1990), then the temperature fluctua- 
tion is given by 



where X{ is the normal vector along the line of sight. Because GM&T predict the 
early formation of objects such as quasars, we assume the universe was reionized 
at high redshift and take this into account by introducing a visibility function, 



which measures the fraction of photons present at redshift z that will reach z = 
without undergoing Compton scattering. If we assume that the electrons which 
scatter each photon see (AT/T) = 0, then we can account for reionization by 
inserting f(z) inside the "Sachs- Wolfe" integral, eq. (1). This means that our 
AT/T results will now depend on h = i/o/100^Mpc _1 and the baryon density 
Qft. We have done most of our calculations for Ml}, = 0.1 which is about the largest 
plausible value. For hQ}, = 0.04, AT/T is only 3% larger on COBE scales. 

We evolve the source fields according to the field equations of motion in the 
Freedman-Robertson- Walker background and calculate the metric perturbations 
due to the energy-momentum of the scalar field by solving the linearized Einstein 
equations. The field equations are 



where v is the vacuum expectation value of the field, and p runs from 1-3 for 
monopoles and 1-4 for textures. Because the defect core size is very much smaller 
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,MMl-(l+ 2 ) 3 / 2 )/21.7 



(2) 




(3) 



3 



than the our grid spacing, we can use the equation for the A — > oo and improve 
the dynamic range of the calculations (Bennett and Rhie, 1990), 



S„-^|(f + ^'-VVU0. (4) 



v 2 J \ a 

The energy-momentum tensor for the scalar fields is given by 

1 'J 1 / -A 2 

e„ ( ,=^ +5 H , 

Ooi =<j) ■ di<j) , (5) 
- 1 / 1 ^ 1 / -a 2 \ 

e y =d i( p ■ + -s^ - - (v0j J , 

and the linearized Einstein equations are (Veeraraghavan and Stebbins, 1990) 

• 2 

h+-h + 3%S C = -8tt (0oo + 6) , (6) 
a cr 

hij - V 2 hij + 2-hij - \didjh + -5ijV 2 h 

_ ° _ 2 - - (7) 

+did k h jk + djd k h ik - -5ijd k dih k i = 167r©y , 

where /i and refer to the traces of /iy and 0y. Thesis used to refer to the 
trace-free terms hij = — \$ijh and 0y = 0y — \&ij®- The perturbation in cold 
dark matter 5 C obeys 

Sc = -\h , (8) 
and the following (non-dynamical) constraint equations must be satisfied, 

1 d i d j hij + \v 2 h = 87r0oo + 3 (-) 6 c --h, 

(9) 



2 ' J lJ 3 V (/ J a 

2 3 v 3 



-djhij — —dih = 87r0Qi 



In theories with "external" sources such as global monopoles or textures, these 
equations impose important constraints on the initial conditions. In particular, 
since ©oo cannot vanish if we are to have interesting density perturbations, the 

— * 

initial energy density fluctuations in the field must be "compensated" by fluctu- 
ations of the opposite sign in the other fields. For the calculations reported in this 
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paper, we have taken h(r) ) = 0, h(r) Q ) = 0, /iy(r? ) = 0, hij(r) Q ) = 0, (f)(r) Q ) = 0, 
and S c (r]o) = — (87r/3)(a/d) 2 9oo(^7o)- Once the constraint equations, (9), are satis- 
fied initially, they will be satisfied at subsequent times if the equations of motion, 
(6)- (8), are satisfied. Extreme care must be taken when evolving these equations 
numerically, because small numerical errors can excite growing mode solutions on 
scales outside the horizon where the physical modes do not grow. When these 
scales finally come inside the horizon, the errors can have grown large enough to 
compete with the physical perturbations that are generated by the source inside 
the horizon. We have found that we can keep these errors under control even when 

— * 

4> takes random values on the scale of 1 grid spacing by going to extremely small 
time-steps (At ps 0.01 Ax). In order to satisfy the second equation in (9), we find 
that global monopole and texture fields must be sufficiently smooth, and this limits 
us to an initial horizon size of ^ 8 grid spacings. 

Another difficulty with evolving eqs. (6)- (8) numerically is that the mixed 
partial derivatives in (7) make it difficult to come up with an explicit differencing 
scheme that is stable. Implicit differencing schemes would be prohibitively expen- 
sive in computer time because we would still have to take very small time-steps 
to avoid the super-horizon scale growing mode solution discussed above since the 
growing mode is a physical mode, not a purely numerical one. 

Fortunately, we have found it possible to remove the offensive mixed partial 
into (7). This substitution yields 



v 

hij =V 2 hij — 2-hij — didjh + 167r<9j0 • dj<j) 
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V=Vo 



+Sij 



V 2 /i+— +4 - S c ---h 
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32n J dr)4> ■ didj(f 
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for eq. (7). The mixed partial derivatives of h and 4> in eq. (10) do not give rise to 
numerical instabilities because the equations of motion for h and <fi ((6) and (4)) 
do not contain any mixed partial derivatives. 

Our numerical simulations use a modified second order leapfrog scheme to 
evolve eqns. (4), (6), (8), and (10) in time. AT/T is determined by integrating 
eq. (1) (with the visibility function (2) inserted) along photon trajectories that 
converge to a point at the end of the simulation. We have done about 25 simulations 
on 64 3 grids and 12 simulations on 100 3 grids for each of the global monopole and 
texture models. The RMS AT/T values for the 100 3 grids are 17.80±2.00Gt4 and 
10.29 ± 1.43Gv 2 for monopoles and textures respectively. (v m and v t refer to the 
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vacuum expectation values of the monopole and texture fields.) These numbers 
can be compared to the RMS fluctuation measured by COBE (AT/T)rms — 
1.10 ± 0.18 x 10~ 5 to yield Gv 2 n = 6.18 ± 1.23 x 10~ 7 for global monopoles or 
Gv\ = 1.07±0.23 x 10~ 6 for textures if we assume that monopoles or textures are 
responsible for the AT/T observed by COBE. 

Figs. 1 and 2 show the simulated full-sky temperature maps for the global 
monopole and texture models respectively smoothed to the same 10° scale as the 
COBE maps. The scale on these plots is given in terms of Gv 2 . Figs. 4 and 5 show 
histograms of the AT/T distributions in the monopole and texture cases. Note 
that after the smoothing the maps contain only about 400 independent pixels. 
Thus, the departures from a Gaussian distribution are not significant. 

Fig. 3 shows the angular power spectrum, AT 2 , of our simulated AT/T maps. 

AT 

— (9,<f ) ) = Y,ai m Y lm (0,<P) , (11) 

l,m 



m 

The solid and dashed curves give the best fit to the predicted power spectrum for 
Harrison-Zel'dovich (n = 1) primordial adiabatic density perturbations (Bond and 
Efstathiou, 1987), 

AT 2_ (n , 2 (2/ + l)r(/ + (n-l)/2)r((9-n)/2) 

1 ~ ^ r ™- ps > 5 r(Z + (5-n)/2)r((3 + n)/2) ' 1 ' 

The best fit of the form, (13), gives n — 1.1 ± 0.3 for global monopoles and 
n = 1.2 ±0.2 for textures when we remove the quadrupole from the fit as was done 
for the COBE data. (The error bars in these fits reflect mainly systematic errors.) 
The Harrison-Zel'dovich value, n — 1, gives a good fit to our simulations and to 
the COBE data, so it makes sense to compare the fit amplitudes. The COBE 
value is Q rms -ps — 6-11 ± 1-46 x 10~ 6 , and we obtain Q rms -ps = 8.7 ± l.GGv^ 
for monopoles and Qrms-PS — 4.7 ± O.hGv 2 for textures. A comparison of these 
values gives Gv^ = 7.0 ± 2.1 x 10~ 7 and Gv 2 = 1.30 ± 0.34 x 10~ 6 consistent with 
the values obtained above. 

It is worth noting that topological defect theories generically predict AT/T 
and Sp/p spectra that are slightly steeper than Harrison-Zel'dovich on very large 
scales (<: 1000 Mpc). The reason for this is that with topological defects, unlike 
inflation, the gravitational fields are generated inside the horizon, so that scales 
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near the horizon have yet to receive their full "share" of perturbations. Thus, the 
largest scales are expected to have less power than the scale-free spectrum (i.e. 
n > 1). (This effect is partially compensated for by the effects of reionization 
which reduce AT/T on small scales.) For cosmic strings, the effect should be even 
more pronounced since the coherence scale of the strings is smaller than that of 
global monopoles or textures (Bennett, Stebbins, and Bouchet, 1992). In contrast, 
power law inflationary models (Davis, et al, 1992) which are able to fit the COBE 
results with a reasonable value for b 8 predict n < 1. Thus, if the four year COBE 
results converge to n — 1.50 ± 0.25, it will be fair to conclude that topological de- 
fect models are preferred over quantum fluctuations during inflation as the source 
of the primordial density perturbations. 





Monopoles 


Textures 


h 




b 8 


Gvf 


b 8 


0.5 


2.49 x 10~ 6 / & 8 


4.03 ±0.80 


4.56 x 10~ 6 / & 8 


4.27 ±0.98 


0.6 


2.19 x 10~ 6 / & 8 


3.54 ±0.70 


3.93 x 10~ 6 / & 8 


3.68 ±0.84 


0.7 


1.96 x 10~ 6 / & 8 


3.17 ±0.63 


3.48 x 10~ 6 /fr8 


3.25 ±0.75 


0.8 


1.82 x 10~ 6 / & 8 


2.94 ±0.58 


3.17 x 10~ 6 / & 8 


2.97 ±0.68 


1.0 


1.64 x 10~ 6 / & 8 


2.65 ±0.53 


2.75 x 10" 6 /^8 


2.57 ±0.59 



Table 1. The scalar field normalizations and best fit bias parameters, b 8 , to the 
COBE-DMR RMS anisotropy at 10° are tabulated as a function of h. la errors 
are reported. 

In order to translate our results into limits on theories of cosmic structure 
formation, we must normalize Gv 2 to give a reasonable spectrum of density per- 
turbations. Because of uncertainties in the relationship between the number den- 
sity of galaxies and the mass density, this normalization is conventionally given in 
terms of a bias factor, b 8 , such that the RMS mass density fluctuation = l/b 8 after 
smoothing with an 8/i _1 Mpc radius top hat. Table 1 gives these normalizations as 
determined in Bennett, Rhie, and Weinberg (1992). 

Table 1 also lists the predicted b 8 values as a function of h = //o/lOO^Mpc- 1 
as determined by a comparisons of the predicted RMS fluctuation at 10°. We 
can see that the predicted b 8 values run very high for small values of the Hubble 
constant. Choosing a large Hubble constant in the Q = 1 universe that we have 
assumed is problematic due to the implied very short age for the universe, but with 
h — 0.7 (the smallest value that is consistent with most of the measured values 
of h) b 8 = 2.5 is within ~ la of the mean in the both the global monopole and 
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texture models. Thus, global monopoles and textures with a bias factor in the 
range of 2.5 — 3 seem to be in reasonably good agreement with the COBE data. 
In a separate study of large scale structure in the global monopole and texture 
models, (Bennett, Rhie, and Weinberg, 1992), we find that these high bias global 
monopole and texture models do reasonably well in matching the observed large 
scale structure. Further work is required to see if the required biasing can be 
obtained dynamically, however. 

Finally, let us compare our results to those of other, well motivated theories 
of large scale structure formation. We find that for similar values of h and bs, 
global monopole and textures predict AT/T on COBE scales that is a factor of 
4-5 larger than the standard CDM prediction (Bond and Efstathiou, 1987). With 
a reasonable bias value, 1.5 ^ b$ 2.5, this model is inconsistent with the COBE 
measurement for h > 0.5 and but perhaps barely consistent for h — 0.5. If we 
demand that 1.5 ^ &8 ~ 2-5 for global monopoles and textures, we find consistency 
with COBE for the entire range, 0.5 < h < 1.0 at the 2a confidence level and 
0.7 < h < 1.0 at the la level. The power law inflationary models discussed by 
Davis, et ai, (1992) can be made consistent with reasonable bg values because they 
have significant contributions to AT/T from tensor modes that do not contribute 
to Sp/p. These models do seem to have some difficulty with forming galaxies early 
enough, however (Adams, et ai, 1992). Other models in with a smaller amount of 
CDM, such as hot + cold DM models or low Q, models seem to fit the COBE results 
reasonably well (Wright, et al, 1992), but they are less attractive theoretically. 
Cosmic Strings + HDM seem to be a good fit to the COBE results (Bennett, 
Stebbins, and Bouchet, 1992), but the theoretical error bars are presently rather 
large. Thus, none of the best motivated models are singled out by the COBE 
results, but global monopoles and textures are arguably the best fit to COBE 
among the £1=1 pure CDM models. Further work into the details of galaxy 
formation and AT/T on smaller angular scales is certainly warranted. 

Acknowledgements: We would like to thank A. Stebbins and D. Weinberg for help- 
ful discussions. This work was supported in part the U.S. Department of Energy at 
the Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48 
and by the NSF grant No. PHY-9109414. 



8 



REFERENCES 



Abbott, L., and Wise, M., 1984, Nucl. Phys. B 244, 541. 

Adams, F. C, Bond, J. R., Freese, K., Frieman, J. A., and Olinto, A. V., 1992, CITA 
preprint. 

Amaldi, U. et al, 1991, Phys. Lett. B 260, 447. 

Bartlett, J., Gooding, A. K., and Spergel, D. N., 1992, Berkeley preprint. 
Bariola, M., and Vilenkin, A., 1989, Phys. Rev. Lett. 63, 341. 

Bennett, D. P., Stebbins, A., and Bouchet, F. R., 1992, IGPP preprint UCRL-JC- 

110803, submitted to Astrophys. J. Lett. 
Bennett, D. P., and Rhie, S. H., 1990, Phys. Rev. Lett. 65, 1709. 
Bennett, D. P., Rhie, S. H., and Weinberg, D. H. 1992, in preparation. 
Bond, J. R., and Efstathiou, 1987, Mon. Not. R. Astron. Soc. 226, 655. 
Cen, R. Y., Ostriker, J. P., Spergel, D. N., and Turok, N., 1991, Astrophys. J., 383, 1. 
Davis, R., Hodges, H., Smoot, G. F., Steinhardt, P. J., and Turner, M. S., 1992, 

preprint. 

Gooding, A. K., Spergel, D. N., and Turok, N., 1991, Astrophys. J., 372, L5. 

Gooding, A. K., et al, 1992, Astrophys. J., 393, 42. 

Park, C, Spergel, D. N., and Turok, N., 1991, Astrophys. J., 372, L53. 

Sachs, K., and Wolfe, A. M., 1967, Astrophys. J., 147, 73. 

Smoot, G., et al, 1992, COBE preprint. 

Turok, N., 1989, Phys. Rev. Lett. 63, 2625. 

Turok, N., and Spergel, D. N., 1990, Phys. Rev. Lett. 64, 2736. 

Veeraraghavan, S., and Stebbins, A., 1990, Astrophys. J., 365, 37. 

Vilenkin, A., 1980, Phys. Rev. Lett. 46, 1169, 1496(E). 

Wright, E., et al, 1992, COBE preprint. 

Zel'dovich, Y. B., 1980, Mon. Not. R. Astron. Soc. 192, 663. 



9 



FIGURE CAPTIONS 



1) A AT/T map for the global monopole model generated by our global field 
evolution code is displayed in a full-sky equal area projection. The scale is 
given in units of Gv^- 

2) A AT/T map for the texture model generated by our global field evolution 
code is displayed in a full-sky equal area projection. The scale is given in 
units of Gvf . 

3) The average AT/T power spectrum is plotted for 12 100 3 monopole simula- 
tions and 12 100 3 texture simulations. The solid and dashed curves give the 
best fit to the power spectrum derived for a Harrison-Zel'dovich spectrum 
of primordial adiabatic density perturbations. The error bars give the RMS 
deviation from the mean, so they reflect the expected deviation for a single 
realization. 

4) A histogram of the pixels in convolved AT/T map shown in Fig. 1. The 
bins on the edges of the histogram include all the points outside the limits of 
the figure. The dashed curve is the histogram for a Gaussian with the same 
RMS. 

5) A histogram of the pixels in convolved AT/T map shown in Fig. 2. The 
bins on the edges of the histogram include all the points outside the limits of 
the figure. The dashed curve is the histogram for a Gaussian with the same 
RMS. 
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